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I Abstract. We present results using a new scaling variable, ^„ in modeling electron- and neutrino- 

nucleon scattering cross sections with effective leading order PDFs. Our model uses all inelastic 
charged lepton F2 data (SLAC/BCDMS/NMC/HERA), and photoproduction data on hydrogen and 
deuterium. We find that our model describes all inelastic scattering charged lepton data, the average 
of JLAB resonance data, and neutrino data at all Q^. This model is currently used by current neutrino 
oscillation experiments in the few GeV region. 
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^ , The field of neutrino oscillation physics has progressed from the discovery of neutrino 

QQ ' oscillation ^ to the era of precision measurements of mass splitting and mixing angles. 

O . Currently, there are only poor measurements of differential cross sections for neutrino 

interactions in the few GeV region. This results in large systematic uncertainties in the 
^ extraction of mass splitting and mixing parameters (e.g. by the MINOS, NOvA , K2K 

p ■ and T2K experiments). Therefore, reliable modeling of neutrino cross sections at low 

' ! energies is essential for precise (next generation) neutrino oscillations experiments. In 

Qj[ the few GeV region, there are three types of neutrino interactions: quasi-elastic, reso- 

^ ■ nance, and inelastic scattering. It is very challenging to disentangle each contribution 

separately, especially, resonance production versus deep inelastic scattering (DIS) con- 
tributions. There are large non-perturbative QCD corrections to the DIS contributions in 
■ this region. 

9^ , Our approach is to relate neutrino interaction processes using a quark-parton model 

to precise charged-lepton scattering data. In a previous communication M\, we showed 
that our effective leading order model using an improved scaling variable describes 
all deep inelastic scattering charged lepton-nucleon scattering data including resonance 
data (SLAC/BCDMS/NMC/HERA/Jlab) [3, 4] from very high to very low (down 
to photo-production region), as well as high energy CCFR neutrino data [5]. 

The proposed scaling variable, is derived using energy momentum conservation, 
assuming massless initial state quarks bound in a proton of mass M. 

, _ 2xiQ^+Mf^+B) 



Q\\ + Vl + (2Mx)2/e2] +2Ax' 



here, is the final quark mass ( zero except for charm-production in neutrino pro- 
cesses). The parameter A accounts for the higher order (dynamic higher twist) QCD 
terms in the form of an enhanced target mass term (the effects of the proton target mass 



are already taken into account using the exact form in the denominator of ). The pa- 
rameter B accounts for the initial state quark transverse momentum and final state quark 
effective AMf^ (originating from multi-gluon emission by quarks). This parameter also 
allows us to describe the data also in the photoproduction limit (all the way down to 

A brief summary of our effective leading order (LO) model is given as follows; 

• The GRV98 LO PDFs M are used to describe the F2 data at high region. 

• The scaling variable x is replaced with the improved scaling variable ^„ (Eq.[l]). 

• All PDFs are modified by K factors to describe low data in the photoproduction 
limit. 

KUQ') = K^alenceiQ') = [1 - GUq')] , (2) 

where Go =1/(1 + 2^/0.71)^ is the proton elastic form factor. At low Q^, [1 — 
Gj){Q^)] is approximately 2^/(2^ + 0.178). Different values of the K factor are 
obtained for u and d quarks 

• The evolution of the GRV98 PDFs is frozen at a value of = 0.80 which is the 
minimum value of this PDFs. Below this Q^, is given by; 

F2{x, < 0.8) = K{Q^) X F2{^ , = 0.8) (3) 

• Finally, using these effective GRV98 LO PDFs {^w) we fit to all inelastic charged 
lepton scattering data (SLAC/BCDMS/NMC/Hl) and photoproduction data on 
hydrogen and deuterium. Note that no resonance data is included in the fit. We 
obtain excellent fits with; A=0.538, 5=0.305, C^i=0.202, Q=0.291, C4=0.255, 
Q=0.189, Cfi=0.621, q'i=0.363, and X^/DOF =1874/1574. Because of the K 
factors to the PDFs, we find that the GRV98 PDFs need to be multiplied by a factor 
of 1.015. 

The measured structure functions data are corrected for the relative normalizations 
between the SLAC, BCDMS, NMC and HI data. The deuterium data are corrected for 
nuclear binding effects [7]. We also add a separate charm pair production contribution to 
lepton-nucleon scattering using the photon-gluon fusion model. This component is not 
necessary at low energies, and is only needed to describe the highest v HERA F2 and 
photoproduction data (since the GRV98 LO PDFs do not include a charm sea). 

Our effective LO model describes various DIS (SLAC, BCDMS, NMC, and HERA) 
and photo-production data down to the = limit. Fig.[l]shows some of comparisons. 
Furthermore, based on duality arguments [Isl], it appears that this model also provides 
a reasonable description of the average value of F2 for SLAC and Jlab data in the 
resonance region, as shown in Fig.|2tleft] . Although no resonance data has been included 
in our fit, our model gives a good description of the most recent 2xFi electron-proton 
Rosenbluth separated data in the resonance region from Jefferson Lab Hall C E94-110 
Collaboration , and also data from the first phase of the JUPITER program at Jlab. Our 
predictions for 2xFi are obtained using our F2 model and Rigg^ [ 10]. For heavy targets, 
nuclear effects are important, especially at low Q^. Recent results from Jlab indicate that 
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FIGURE 1. Comparisons of the predictions of our effective LO model for F2 to charged lepton inelastic 
scattering data, [left] DIS F2 proton data (SLAC, BCDMS, NMC), [right] HI F2 proton. 
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FIGURE 2. Comparisons of the predictions of our effective LO model to resonance electro-production 
data on protons (which was not included in our fit). Shown are F2 proton [left], and 2xFi proton data [right] 
(from the Jefferson Lab Hall C E94-1 10 Collaboration). The predictions for 2x^1 are obtained from our 
model for F2 with Rigga- In the right plot, the solid line uses GRV98 PDFs, and the dashed line is our 
previous model using GRV94 LO PDFs. 



the Fe/D ratio in the resonance region is the same as the Fe/D ratio from DIS data for 
the same value of ^ (or ^w)- Future Jlab experiments with deuterium and heavy nuclear 
targets (e.g. JUPITER) will provide a high statistics data in the resonance region which 
will be very important to improve our model at region. 

In neutrino scattering, in addition to the vector structure function, there is an axial 
vector structure function contribution. At the = limit, the vector structure function 
goes to zero, while the axial-vector part has a finite contribution. At high Q^, these 
two structure functions are expected to be same. Thus, it is important to understand the 
axial- vector contribution at low by comparing to future low energy neutrino data (e.g. 
MINERvA [14]). As a preliminary step, we compare the CCFR and CDHSW [11] high 
energy neutrino data with our model, assuming that the vector contribution is the same 
as the axial vector contribution. We find that the CCFR/CDHSW neutrino data are well 
described by our model. 

We are currently working on constraining the low axial vector contribution using 
low energy CDHSW and CHORUS JH data. The form of the fits we plan to use is 



motivated by the Adler sum rule Ill3ll for the axial vector contribution as follows: 



K.ea-aAQ') = €^^: K.alenceiQ^) = [I - FM)] ( . (4) 
+ Cls-ax \Q + Clv-flx / 

where Fa{Q^) = —1.267/(1 + 2^/1.00)^. Nuclear effects for heavy target are also 
important and may be different for the vector and axial vector structure functions. Future 
measurements on the axial vector contribution from the MINERvA experiment lfl4l] will 
be important in constraining this model. 
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